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Rashba spin-orbit coupled semiconductor-superconductor hybrid structures in the presence of 
Zeeman splitting have emerged as the first experimentally realizable topological superconductor 
supporting zero-energy Majorana bound states. However, recent experimental studies in these hy- 
brid structures are not in complete agreement with the theoretical predictions, for example, the 
observed height of the zero-bias conductance peak (ZBCP) associated with the Majorana bound 
states is less than 10% of the predicted quantized value 2e^//i. We try to understand the sources of 
various discrepancies between the recent experiments and the earlier theories by starting from a mi- 
croscopic theory and studying non-equilibrium transport in these systems at arbitrary temperatures 
and applied bias voltages. Our approach involves quantum Langevin equations and non-equilibrium 
Green's functions. Here we are able to model the tunnel coupling between the one-dimensional 
semiconductor-superconductor hybrid structure and the metallic leads realistically; study the role 
of tunnel coupling on the height of the ZBCP and the subgap conductance; predict the nature of the 
splitting of the ZBCP with an increasing magnetic field beyond the critical field; show the behavior 
of the ZBCP with an increasing gate-controlled onsite potential; and study the evolution of the 
full differential conductance across the topological quantum phase transition. When the applied 
magnetic field is quite large compared to the Rashba splitting and the bulk energy gap is much 
reduced, we find the ZBCP even for an onsite potential much larger than the applied magnetic field. 
The height of the corresponding ZBCP depends on the tunnel coupling even at zero temperature 
and can be much smaller than 2e'^ /h. 



I. INTRODUCTION 

Low-dimensional conventional-superconductor setups 
have attracted growing attention recently, particularly 
in the context of the study of fluctuations^ and quantum 
phase transitions.^ On the other hand, Kitaev^l was the 
first to propose the idea of realizing Majorana fermions 
as localized states at the ends of an ideal ID (spinless) 
p-wave superconducting wire motivated by the search for 
robust quantum information storage. This, eventually, 
generated a growing interest to produce emergent Majo- 
rana fermions in the solid-state laboratory setting.l3HIl A 
particularly promising proposal to engineer effective p- 
wave superconductivity is by using a Rashba spin-orbit 
coupled semiconductor in proximity with a conventional 
s-wave (spin-singlet) superconductor and in the presence 
of Zeeman splitting due to an applied magnetic field.'^HiSl 
The Rashba coupling creates two helical bands which 
wind counterclockwise and clockwise. Introduction of the 
proximity induced superconductivity gives rise p±ip pair- 
ing for the two helical bands. One then needs to apply a 
magnetic field to break the time-reversal doubling of the 
fermionic states on the spin-orbit coupled semiconductor. 
This thus leads to a topologically nontrivial single species 
p+ip superconductor as originally conceived by Kitaev.l^l 
In the presence of a magnetic field, there is an effective 
s-wave pairing between the two helical bands apart from 
the p±ip pairing for each band. The strength of the p- 
wave component of pairing depends on the amount of 
Rashba splitting. Parameters are to be chosen so as to 



keep a finite p-wave component of the induced pairing 
gap such that the Majorana end states in a wire geome- 
try are protected from the bulk excitations. 

The recent experimental efforts to detect Majo- 
rana quasiparticles in one-dimensional semiconductor- 
superconductor hybrid structures in the presence of 
strong spin-orbit coupling and a Zeeman fielcP^HUl dis- 
play sever al disc repancies with the existing theoretical 
predictions.'^i^ These experimental studies have mostly 
concentrated in measuring a zero-bias conductance peak 
in the tunneling differential conductance of hybrid struc- 
tures at low temperatures by tuning the applied mag- 
netic field along the one-dimensional (ID) structure. It 
has been theoretically predicted that the existence of 
zero-energy Majorana bound states (MBSs) in these sys- 
tems manifests through a quantized zero-bias conduc- 
tance peak (ZBCP) of height 2e^//i at zero-temperature 
and above a critical magnetic field {3^) which is de- 
termined by the gating of the spin-orbit coupled semi- 
conductor and the magnitude of the induced supercon- 
ducting gap. On the other hand, the height of the 
observed ZBCP in the experiments is much lower (less 
than 10%) than the theoretical predictions. It has been 
also seen that the measured ZBCP splits at higher mag- 
netic fields beyond Be- There are also additional dis- 
crepancies between theory and experiments which are 
not yet completely resolved. For example, it was the- 
oretically predicted that the superconducting gap would 
close across a magnetic-field or onsite-potential driven 
topological quantum phase transition between a non- 
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topological (s-wave) superconducting phase and a topo- 
logical (p-wave) su percon ducting phase with zero-energy 
Majorana modes P *^^ l ^^ l However, some experiments^^ 
have claimed to observe such gap closing and some have 
not seen a gap closingli^ The expected disappearance of 
the ZBCP as the Zeeman field is rotated from the wire 
axis in the plane formed by the wire axis and the di- 
rection of the effective spin-orbit field has not been ob- 
served in the latest e xperim ent!^ but it was seen in the 
previous experiments.l^IlSI The emergence of real MBSs 
in these systems can also be probed by looking for a 
fractional Josephson effect in some type of interference 
experiments,-—"' which manifests though a Att period- 
icity in an ac Josephson measurement. These type of 
doubled-Shapiro-step measurementJ^ll are complex and 
harder to interpret theoretically with a realistic micro- 
scopic modeling. 

The existing theoretical studies to calculate 
current-voltage characteristics of various topological- 
superconductor configurations use tools such as 
formal isms bas ed on the Landauer-Biittiker scattering 
theory,'Sl22l2IIIIl or the Keldy sh no nequilibrium Green's 
function (NEGF) formalisir|2Sl2l (as well as other 
variant approaches, see for instance Ref. |27j). With the 
latter, models have been more idealized, while with 
the former, models do not take into account the bath 
or the coupling with it explicitly but in an effective 
manner. While a bath can induce decoherence in the 
semiconductor-superconductor hybrid structures, (and 
thus limit the performance as a qubit of a pair of MBSs) , 
its coupling to the hybrid structures also controls the 
height of the ZBCP at finite temperature. Therefore, 
an explicit modeling of the baths and their couplings 
is an important issue to understand the recent experi- 
mental data. Recently, we have extended the quantum 
Langevin equations and Green's functions (LEGF) 
formalism to study nonequilibrium transport in p-wave 
superconductors.^^ One advantage of the LEGF method 
is that it starts with an explicit Hamiltonian for the bath 
and the coupling along with the system Hamiltonian, 
thus the role of the bath and the coupling is quite clearly 
elucidated. In this paper we provide a detail derivation 
of the LEGF method and further extend it for the case 
of spin-orbit coupled semiconductor-superconductor 
heterostructures in the presence of an external magnetic 
field. We demonstrate a nonmonotonic dependence of 
the height of the ZBCP with the applied magnetic field 
and the onsite potential. We also show the oscillation 
of the splitting width in the ZBCP with an increasing 
magnetic field beyond Be- We further show how the 
increase of the tunnel coupling can modify the height 
of the ZBCP at finite temperatures and control the 
appearance of sub-gap conductance. In short, we provide 
a complete theory of linear and nonlinear transport in 
the topological superconductors at all temperatures and 
elucidate many interesting results for the dependence of 
the ZBCP on various experimentally tunable parameters. 
The rest of this paper is organized as follows: in Sec. In] 



we introduce the model of the wire, in Sec. Ill we explain 



the treatment using quantum Langevin equations and 
Green's functions, and in Sec. |IV]we present the results 
mentioned above. 



II. MODEL AND MAJORANA MODES 

Let us consider a single-channel semiconductor 
nanowire with a strong Rashba spin-orbit coupling (for 
example, as in InSb or InAs) in close proximity to an 
ordinary (s-wave) superconductor (such as, for instance, 
NbN, NbTiN, or Al) and in the presence of an external 
magnetic field applied along the axis o f the n anowire. 
It has been theoretically predicted earlieil^Ei^ that this 
system undergoes a topological quantum phase transi- 
tion at a certain critical magnetic field, Be = \/ + 
where A is a proximity induced superconducting gap and 
/i is an onsite potential for the nanowire. For an applied 
magnetic field B > B^ the hybrid structure is driven 
into a chiral p-wave topological superconducting phase 
supporting two real zero-energy MBSs at the two ends 
of the nanowire. The requirement of finite p-wave com- 
ponent of induced pairing gap for the protection of the 
MBSs along with the relation for the critical field B^ 
impose a stringent restriction on the onsite potential or 
carrier density. One needs a low carrier density in or- 
der to satisfy those conditions, therefore the engineered 
structures are susceptible to disorder. The ZBCP in the 
differential conductance of the tunneling transport cal- 
culations is a signature of the emergence of MBSs in the 
topological superconducting phase. The semiconducting 
wire in the semiconductor-superconductor heterostruc- 
ture is modeled by a tight-binding lattice Hamiltonian 
of electrons with a prox imity induced s-wave BCS pair- 
ing of amplitude A.EnH The wire has N lattice sites and 
each of its two end sites is coupled to an infinite metallic 
bath which is itself modeled by a one-dimcnsional tight- 
binding system of free electrons. The effective Hamilto- 
nian of the system consisting of the nanowire ,1^'^° the 
baths, and the tunnel couplings reads 
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Here a^^. and a™J, denote respectively an electron cre- 
ation operator on the semiconductor nanowire and on 
the m*'' bath (here m = L,R). The Hamiltonian of the 
nanowire is denoted by Hw, that of the to"^ bath by "Hg , 
and the tunnel coupling between the nanowire and the 
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The hopping amplitude of an electron 
on the wire is 7, fj, is an onsite potential energy of the 
wire controlled by gating, a is the strength of the Rashba 
spin-orbit coupling on the semiconductor wire and B is 
an applied magnetic field along the axis of the wire. We 



shall choose below the values of the different parameters 
for our plots to be consistent with the values reported in 
recent experiments. 

Next we introduce the following local transformation 
for the electron operator of the semiconductor nanowire, 
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where c^,/,^ is a Majorana (or real) fermion opera- 
tor. The semiconductor-superconductor hybrid struc- 
ture Hamiltonian in the Majorana- fermion basis is thus 
rewritten as: 
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A. Low Energy Spectrum 

When B > Be we find that the above 47V x AN ma- 
trix of the Hamiltonian in Eq. ([3| has two degenerate 
zero-energy eigenstates for large N. These zero-energy 
Majorana modes are separated from all the other energy 
eigenvalues by an energy gap of order A. However, there 
is a splitting between the two Majorana modes for smaller 
values of N and the energy of the Majorana modes is not 
exactly zero. The crossover between these two, longer- 
vs smaller-length, behaviors of the MBSs is determined 
by the coherence length (which is roughly given by 
ttvf/A, where vp is the Fermi velocity) of the spin-orbit 
coupled semiconductor-superconductor hybrid structure. 
The coherence length is longer for a smaller proximity- 
induced superconducting gap or spin-orbit coupling, and 
it also depends on the applied magnetic field and the 
onsite potential.'^ It is believed that the length of the 
nanowires in recent experi ment s was of the same order 
as the coherence length ^gj ^^ l ^il Therefore it is expected 
that there is an overlap between the two MBSs at the two 
ends of the nanowire in these experiments. The overlap 
between the MBSs at the two ends goes to zero for a long 
nanowire where the length of the wire is much longer than 
the coherence length ^o- Whenever there is splitting in 
the energy of the MBSs of an isolated hybrid structure, 
it is expected to show up as a splitting in the ZBCP of 
tunneling measurements as long as the broadening of the 
conductance peak due to the coupling with the bath and 
the temperature is smaller than the energy splitting of 
the two MBSs. We also find that as the applied mag- 
netic field is increased much above B^, the bulk energy 
gap is substantially reduced. There are still zero-energy 



real eigenstates at these parameters, however these zero- 
energy states are very delocalized and extend into the 
bulk of the wire. When the onsite potential is increased 
at a higher magnetic field, we find zero-energy states in 
the isolated Hamiltonian of Eq. ([3| even when the ap- 
plied field becomes smaller compared to the critical field 
at that value of the potential (cf.'22). But when the onsite 
potential is substantially increased the spectrum of the 
Hamiltonian in Eq. ([s]) becomes again gaped near zero 
energy, and there is no zero-energy state (same as for the 
simpler Kitaev model). Below, we consider this regime 
of parameter sets carefully in our analysis of transport, 
because we suspect that the recent experiments might 
be probing it. Although the spectrum and local density 
of states of this regime had already been analyzed via 
exact diagonalization of the isolated-wire Hamiltonian, a 
transport analysis had not been carried out previously. 



III. QUANTUM LANGEVIN EQUATIONS AND 
GREEN'S FUNCTION FORMALISM 

Here we develop a steady-state non-equilibrium trans- 
port theory for the spin-orbit coupled semiconductor- 
superconductor hybrid structures by employing quantum 
Langevin equations and Green's functions. We obtain a 
set of generalized quantum Langevin equations of mo- 
tion for the nanowire's operators in the Majorana basis 
following Refs..'22H2Sl The description here for the hybrid 
structures is similar to that for the Kitaev chain in our 
earlier study!^- We assume that the metallic baths at the 
two ends of the nanowire are disconnected from the wire 
for all times t < Iq. Each bath is in thermal equilibrium 
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at a specified temperature and a chemical potential 
/im for both spin components of the bath. We connect 
the baths to the nanowire at a time s-nd we are inter- 
ested in the steady-state properties of the nanowire that 



set in after the end of any transient dynamics. We write 
the Heisenberg equations of motion for the operators on 
the nanowire in the Majorana basis at i > tg, 
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for / = 1, 2, 3, and with the convention that CAfi.a — quantum Langevin equations (see Appendix|A]for the full 
CA,Af-i-i,cr = cb,o,<t = CB,Ar+i,cr = 0. The Heisenberg equa- 
tions of motion for the bath operators (where a =t, 4, and 
m ^ L, R) are: 



expressions) where we identify 77^,0- as a noise contribu- 
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The equations of motion of the wire operators, Eqs 
in the Majorana basis involve the bath variables a" 
with m = L, R, that we can eliminate by replacing with 
their exact solutions. For that we note that the equations 
of motion of the each bath, given by Eqs. ( 6|7|8 1, are a 
set of linear coupled equations with an inhomogeneous 
part given by 17^0™, o-/^- We solve these equations of 
motion using the single-particle retarded Green's func- 
tion of the isolated baths, which is given by g^^{t) — 
— i0(t)e^*^"*/'' where Hjp is the single-particle Hamilto- 
nian of the spin cr component of the m^^ bath. As here 
the single particle retarded Green's function of the each 
bath is the same for both spin components, hereafter we 
avoid the spin index in 17™+ (i), and write it as g"^^{t). 
One finally finds that the solution for the boundary site 
on the m"^ bath is given by (for t > to) 
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dt' gt+{t-t')'^-^aNAt')- (10) 



Plugging these solutions into Eqs. (4[5) for the wire op- 
erators in the Majorana basis, we get a set of generalized 



tion from the spin-cr component of the m}^ bath while the 
terms involving ^^(i) are the corresponding dissipative 
terms. 
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S+(t) = (^)^9ri+(t), and E,- (t) = [I]+(<)]t.(12) 

The noise depends on the bath's initial distribution 
which we have chosen to correspond to thermal equi- 
librium. The properties of the noise are written most 
conveniently in the frequency domain. Let us consider 
the limit tg ^ —00, and introduce the Fourier trans- 
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{iJf^fgtn^) (here g^+i^) is the m''^ 
bath single-particle Green's function evaluated at the 
first site). We also have cj^ ^ ^(w) = Cfj^„i.(j{—^)- We use 
the definition r,„(a;) = -/toP+J/tt = {i^J pjn{u}) , 
where Pm{^) is the local density of states of either spin 
component at the first site (n = 1) on the m'^ bath. With 
these definitions it is easy to show that the noise-noise 
correlations are given by 

{Vl A^)f}m.a'{'^')) = r;(a;)/(tJ,/l;,Tj)(5(w - Uj')5lm5a,a' , 

(13) 

where /(w,/i„,r,„) = l/{exp[(;iw - fJ-,n)/kBT.m] + 1} 
is the Fermi distribution function and fc^ is the Boltz- 
mann constant. The eigenvalues and eigenfunctions of 
the m*"^ bath Hamiltonian are given by e™ = —27^ cos q 

and ipq{p) — %/2sin((7p), where q lies in the range [0, tt] 
and p is an integer. By convention, the wave functions 
are normalized so that {ipqlipq') = TrS{q — q'). Therefore, 
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After converting the q sum to an integral in the range 
[0,7r], we find that, within the band-width of the bath 
{\huj\ < 27,„), 5"+(a;) is given by 
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Using these single-particle retarded Green's functions of 
the baths, we find 
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We solv e the ge neralized quantum Langevin equations 
in Eqs. ( Al|A2 ) by Fourier transform and we get the 
following steady-state solution for the operators on the 
nanowire. 
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ci{uj) = Gl^{uj)h„i{uj), where 

m — 1 

G+{ij) = Z-\lj), Zi„M = ^v,n{uj) + AiM, (16) 
and we use the following notation for simplicity, ca.i,1; = 

C4/-3, <^B,l,^ = C4;_i, CA,l,i = C4;_2, CB,l,i = C4i for 

I = 1, 2, 3, iV. The expressions for <I>;m(cj), Aim((u;) 
and hjn{uj) are given in Appendix |b[ Here G+(a;) is 
the Green's function of the full system consisting of the 
nanowire and the baths. We can calculate noncquilib- 
rium steady-state properties of the hybrid structures us- 
ing the above solutions for the operators. We now pro- 
ceed in the next section to evaluate the electrical current 
in the system under an arbitrarily large applied bias volt- 
age and finite temperature. 



IV. CURRENT- VOLTAGE CHARACTERISTICS 
AND ZERO-BIAS CONDUCTANCE PEAK 

We define a local charge density at the boundary site 
of the nanowire and find an expression for the electri- 
cal current through the nanowire using the local density 
and the continuity equations. '^'^ ■^^ We call jm{t) the in- 
ward electrical current flowing from the to"^ bath into 
the nanowire. It is given by 
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wher e (. ..) denotes averaging over noise using the result of 
Eq. ( [T3| . Notice that, due to the proximity-induced su- 
perconductivity in the nanowire, the total electron charge 



is not conserved. Thus, the electrical current from the left 
bath into the nanowire is not necessarily equal in mag- 
nitude to that from the nanowire into the right bath, for 
arbitrary chemical potentials in the baths. In order to be 
specific, from now on we explicitly discuss the electrical 
current from the left bath. We find from Eq. [IT) after 
using Eq. [9] with — )■ —00, 

jL(t)=2Im((^ {a\jt)[riLAt) 



+ dt'j:t{t-t')aUt')]})). (18) 

Let us derive each part of the above expression separately. 
After averaging over the noise using Eq. |13[ we find for 
the first part of Eq. [18) 

{J2 «U(i)^L,.(i)) - r ^{[~Gt^i-u^)~^GU~u:) 

+iG+,{-oj) - G3+3(-c.)]ri(c.)/(^,ML,Ti) + hG+ (-^) 
-iGA^uj)+tGU-Lo) - GU-u:)]TLiuJ)fiuJ,^IL,TL)}. 

After applying the local Majorana basis transformation 
we find for the second part of Eq. [18) 
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in where we can now use the solutions of the operators 
on the nanowire given by Eq. 16 For example, the first 
term is given by 

- / di'E+(t-i')( E cAsAtyAsAt')) 
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The noise average in Eq. [19] is carried out using Eq. |B3| 
along with the noise correlation properties from Eq. |13| 
In the steady state, jL(i) is independent of time. Here we 
evaluate the current-voltage characteristics of the hybrid 
nanowire structures in the steady state. We can calculate 



the full Green's function in Eq. 16 numerically and use 
them to find the current at zero or finite temperatures 
using the above results. 



A. Numerical Results 

We now focus on the 'symmetric' case for which the 
two baths are identical (7^ = 7^) and are connected to 
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FIG. 1. Zero temperature dl/dV vs. V with increasing mag- 
netic field. Everywhere A'' — 40, 7m — "/ — 1, a = 0.2, /i = 
and |A| = 0.3. The magnetic field is chosen as follows: (a) 
B = 0.2, (b) B = 0.3, (c,d) B = 0.4, (e) B = 1.25, (f) 
B = 1.5, (g) B = 1.9, and (h) B = 2.1. Everywhere 7^ = 0.2 
(m = L,R), except 7^ = 0.5 in (d). 



the nanowire by identical contacts (7^ = 7^ or Tl — 
F/j = r) while fiL = — —fJ'R and Tl = T = T^. In this 
case the steady-state currents jL,R s^re equal up to a sign 
and 



duj 
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Here T{uj), which contains contributions to transport 
coming from both electrons and holes, can be simply in- 
terpreted in the symmetric case via the zero-temperature 
differential conductance (with /I = eV, e = 1) 



dV djl 
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One can also calculate the finite temperature differ- 
ential conductance from Eq. [20] which is important to 
understand the recent experiments. At zero tempera- 
ture and lower magnetic field (for relatively low /x) the 
height of the ZBCP in the topologically nontrivial phase 
is independent of the contacts with the wire as shown 



in Fig. [T] panels (c,d) for two different contact strengths. 
Those two plots show that the width of the ZBCP in- 
creases with increasing strength of contacts and the con- 
ductance within the pairing gap becomes also finite for 
stronger contacts. The latter can be understood because 
the dephasing induced by the baths becomes substan- 
tial with an increasing strength of the contacts and that 
creates finite conductance even within the superconduct- 
ing pairing gap. However, at finite temperatures such 
that ksT > F, the height of the ZBCP depends on 
the strength of the contacts.!^ When that happens, the 
height of ZBCP falls rapidly with increasing tempera- 
ture of the baths as illustrated in Fig. [2] panels (a,b,d). 
If F > ksT the effect of temperature on the height of 
the peak is minimal and the ZBCP remains almost like 
at zero temperature [see Fig.[2](c)]. 




FIG. 2. Finite-temperature dl/dV vs. V at different temper- 
atures (ksT) and strengths of contacts. Everywhere = 40, 
7m = 7 = 1 (with m ^ L,R), n = 0, a = 0.2, B = 0.4 
and |A| = 0.3. The rest of the parameters are as follows; 
(a) 'y'^ = 0.2, ksT = 0.01; (b) 'y'^ = 0.1, ksT = 0.01; (c) 
7m = 0.5, kgT = 0.01; and (d) 7^ = 0.2, ksT = 0.05. 

In Fig. [1] we show the nature of the zero-temperature 
differential conductance with an increasing magnetic 
field. When the applied field is B < Be (with Be ^ A 
since ^ = in the case of the figure), there is a gap 
in the zero-temperature differential conductance around 
zero voltage for relatively weaker contacts; as is shown in 
Fig. [TJa). The gap in dl/dV closes at B = Be as shown 
in Fig. TJb). For B just above Be, the gap reopens and 
a ZBCP appears in the dl/dV characteristics with its 
height being 2e^ /h at zero temperature for the ideal case; 
this is shown in Fig. [T](c). As we further increase B, the 
pairing gap separating the MBSs from the higher excita- 
tions is reduced. Here we find a splitting in the ZBCP as 
well as finite voltage conductance peaks inside the super- 
conducting pairing gap. The height of the split MBS peak 



is e /h and it is shown in Fig. [l|e). The splitting of the 
ZBCP disappears as we increase B further, however the 
height of ZBCP does not retrieve to its full value of 2e^ /h 



as is shown in Fig. [ijf). As the field is increased even 
more, we find first a new closing of the gap, and finally a 
topologically trivial gapped superconducting phase with- 
out MBSs. These are shown in Fig.[TJg,h). Recently, an 
oscillation in the splitting of the ZBCP with increasing 
magnetic field has been proposed as a smoking gun for 
detection of the MBSs in these hybrid structures!^ Here 
we find from a full microscopic transport calculation that 
one can expect one such oscillation in the splitting of the 
ZBCP with increasing field. However, we also show that 
the height of the ZBCP is reduced with increasing field 
which is an important piece of information for the ex- 
periments, and has indeed been observed recently An 
additional important finding is the retrieval of the gap- 
closing phenomena at higher magnetic fields coincident 
with the disappearance of the MBSs, which is analogous 
to the lower-field scenario just before the emergence of 
the MBSs and reminiscent of the way the topological- 
nontopological quantum phase transition takes place as 
a function of onsitc potential for a Kitacv chain.'SSl 
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FIG. 3. Zero-temperature dl/dV vs. V with decreasing 
Rashba spin-orbit coupling. Everywhere A'^ — 40, 7™, = 
7 = 1, 7^, = 0.2 (with m = L,R), n = 0, B = 0.4 and 
|A| = 0.3. The Rashba spin-orbit coupling is chosen as fol- 
lows: (a) a = 0.1, and (b) a = 0.025. 

Next we study the role of the Rashba spin-orbit cou- 
pling a to observe the ZBCP associated with the zero- 
energy Majorana fermions. Usually these semiconductor- 
superconductor hybrid structures exhibit small Rashba 
splittings which is of order akp ^ 0.1 meV where kp is 
the Fermi momentum for /i = 0. Therefore we look for 
the effect on the ZBCP of a decreasing strength of a. The 
pairing gap around the ZBCP is reduced for smaller a as 
shown in Fig. [sf^a). Further, the gap is fully closed for 
an even smaller value of the spin-orbit coupling and the 
ZBCP associated with the MBSs disappears; as shown in 
Fig-ib). 

The onsite potential of these hybrid structures can be 
tuned by applying gate voltages. However it is difficult 
to measure this potential in the current experiments (si- 
multaneously, the number of channels in the nanowire 
is not exactly known). It has been predicted earlier that 
there can be ZBCP in the confined hybrid structures even 
when n > B for a relatively high magnetic field.'^ Here 
we check the zero-temperature behavior of the ZBCP at 
relatively lower and higher magnetic fields. At a lower 
fixed magnetic field B, there is a ZBCP when /i is such 




FIG. 4. Zero-temperature dl/dV vs. V with increasing /i at 
relatively lower values of B. Everywhere = 40, 7m =7 = 1, 
7m = 0.2 (with m ^ L,R), a = 0.2, B = 0.4 and [Aj = 0.3. 
The gate-controlled onsite potential is as follows: (a) /i = 
0.25, (Be = 0.39); (b) fi = 0.3, (B^ = 0.42); (c) = 0.5, 
{Be = 0.58); and (d) ^ = 1, {Be = 1.04). 



that B > Be = + /i2 as shown in Fig. [ij^a). Only 

one sub-band of the wire is occupied at this value of fi 
and the spectrum of the isolated wire is gapless. For a 
higher value of fi both the sub-bands are occupied, and 
then there are two MBSs coupled by s-wave pairing at 
each end of the wire. The energy of the two MBSs at the 
same end of the wire in the topologically trivial phase 
(B becomes smaller than Be for the corresponding /i) is 
non-zero. Therefore the differential conductance shows 
two finite-voltage strong peaks at these values of ji. All 
these features are shown in Fig. |4|^b,c,d). Interestingly, 
the height of the split peaks remains 2e^ /h which is very 
different from the behavior of the split peaks with in- 
creasing mag netic field at n = [cf. Fig. [TJe)]. While 
the split peaks with height 2e^//i are coming from the 
two coupled MBSs at the same end of the wire, the split 
peaks at /i = with increasing B are due to the overlap 
of two MBSs from the two ends of the nanowire in the 
topologically nontrivial phase. 

On the other hand, the behavior of the ZBCP at zero 
temperature with increasing /i is quite different at a rel- 
atively higher magnetic field. We start with values of fi 
and B so that B > i?c at that /x. There the ZBCP is 
split and the height of the peaks is almost e'^/h as shown 
in Fig. [5ja). As we increase ^ the splitting of the ZBCP 
first disappears and then reappears. These are shown 
in Fig. [sjjbjC). The values of A, fi and B in Fig. [5jb,c) 
correspond to a topologically nontrivial phase. When fi 
reaches a value where B becomes exactly equal to Be 
of the corresponding /i, the superconducting pairing gap 
closes [see Fig. [5|d)]. For an even larger fi the hybrid 
structure would be in principle in the topologically triv- 
ial phase for that B as now B < B^ dX that fx. However, 
we find that the ZBCP reappears again with an increas- 
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(a) (b) 




FIG. 5. Zero temperature dl /dV vs. V with increasing fj, at 
larger B. Everywhere A'" = 40, 7m =7 = 1, 7m = 0.2 (m = 
L,R), a = 0.2, B = 1 and jA] = 0.3. The gate-controlled 
onsite potential is as follows: (a) = 0.25, {Be = 0.39); 
(b) = 0.5, (Be = 0.58); (c) = 0.95, {B^ = 1.00); (d) 
M = 1, [Be = 1.04); (e) = 1.5, (B^ = 1.53); (f) a* = 2.5, 
(Be = 2.52); (g) /i = 2.9, (B^ = 2.92); and (h) ^ = 3.0, 
(Be = 3.01). 



ing but the height is smaller than 2e^//i, which is shown 
in Fig.jsje). As we further increase the onsite potential, 
the ZBCP splits and finally disappears with another gap 
closing step. For a very large value of ^ the spectrum is 
fully gapped at the given field. It is interesting to com- 
pare the emergence of ZBCP in the topologically trivial 
phase with the findings of RefOH for single and multi- 
channel Majorana wires. 

We have also found that the height of the ZBCP at 
high magnetic fields, both in the topologically nontrivial 
and trivial phases, is susceptible to the strength of tunnel 
contacts even at zero temperature. The height increases 
for a stronger tunnel contact with the baths (lower or 
smoother barrier potential) even at zero temperature. 
This behavior is different from the nature of ZBCP in the 
topologically nontrivial phase at a relatively smaller mag- 
netic field [as shown in Fig. [Tjpanels (c,d)]. The splitting 
oscillation of the ZBCP with changing gate-controlled on- 
site energy has been observed in Ref. [TB]where the height 



of the ZBCP is much smaller than 2e^//i even at the low- 
est temperatures accessible to experiments. Therefore, 
we suspect that the recent experiments have probed this 
regime at high B and /i. However, the p-wave compo- 
nent of the induced superconducting state is very small 
for these parameters and the Majorana end states are not 
well protected from the bulk excitations due to a reduced 
p-wave component of the induced pairing gap. 



V. DISCUSSION AND PROSPECTS 

In summary, we have provided a complete theory of 
both linear and nonlinear transport in realistic micro- 
scopic models of hybrid semiconductor-superconductor 
heterostructures with spin-orbit coupling and have shown 
the presence of topological superconducting states as re- 
vealed by tunneling transport. This is manifest in the 
current-voltage characteristics via a ZBCP with a dis- 
tinctive variation as a function of temperature and var- 
ious experimentally tunable parameters. We used and 
provided a detail derivation of the LEGF method, which 
we further extended from our previous worlP^ to include 
the elements of spin-orbit coupling and applied magnetic 
fields. A nonmonotonic dependence of the height of the 
ZBCP with the applied magnetic field and the onsite po- 
tential constitutes a characteristic feature of these sys- 
tems that cannot be explained away by the presence of 
spurious conduction channels (due to the presence of dis- 
order or to other features of the nanofabrication of the 
heterostructures). Moreover, the predicted alternating 
splitting and width modulation in the ZBCP with an in- 
creasing magnetic field beyond is a feature that can- 
not be generically obtained via a two-level system mecha- 
nism (cf. Ref. "^S) , not even with a Zeeman-field tuning of 
single-particle localized levels into degeneracy (similar to 
the scenario for the so called singlet-triplet Kondo effect 
in quantum dots). The above, combined with the char- 
acteristic feature of the coordinated opening and closing 
of the superconducting gap as the MBS's ZBCP first ap- 
pears (and then again when it disappears), constitutes a 
reasonably individualized scenario that would present a 
compelling case for the observation of Majorana states. 
We have further shown how the increase of the tunnel 
coupling acts to modify the height of the ZBCP at fi- 
nite temperatures and controls the appearance of sub-gap 
conductance. Our modeling of the contacts is more realis- 
tic than in more simplified models prevalent in the litera- 
ture and as a result we can meaningfully compare features 
arising due the the presence of MBSs versus other quasi- 
particle excitations. How these characteristic features 
are reflected in other types of transport experiments in- 
volving interferometry in Josephson-j unction geometries 
is an important question for further theoretical and ex- 
perimental study. 
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Appendix A: Generalized Quantum Langevin Equations 

The generalized quantum Langevin equations for the wire operators in the Majorana basis are explicitly given by 
the following two expressions: 

2(m-7) , 2(5 + A) ^ 2(i?-A) 7 7 a 

CA,l,cr — V CB,l.a H V Oa,^CB,l,l H r I C_B / + — ^CB.l + l.a — -^CB,l-l.cr + T"'^,i^B,l+l,f 

n n n n n n 

dt'Y.j^{t - t')-{cA,l,ait') + iCB,l,a{t')) 



iTll^ + i I dt'^l{t - t')\{cA,iAi') - icB,iAt'))) + Si,n{ ~ im,a - i I dt'^+{t - t') 

J to ^ J to 

1 1 

y-^{cA.N,a(t')^iCB^N,a(t'))+^f)\,a+^ \ dt' i:]^{t ^ t') - {cAM ,a{t') ^ iCB M ,a (t'))) , (Al) 



J to 

2(a^-7) 2(i?-A) ^ 2CB + A) 7 7 a 

CB,l,cr — ^ Ca^Im ^ Oa.,t(^A.,l,l ^ <J<7,ICA,L^ + jCA,l + l^a + ^CA.l-l.a — -^0(T,^C^,i + l,t 

Q/ Q/ POO 

~ -^3a,tCA,l-l,i + -^Scr,lCA,l-l,^ + -^S„,^CA,l+l^l + VL,a - J dt'Y.+ {t - t')-{cA,l,a{t') + iCB,!^^')) 

-via- dt'j:i{t-t')-icA,lAt')~iCB,lAt'))) +Sl.N{-VR,a~ / dt'^+{t-t') 

Jto ^ Jto 

1 1 

^^{cA.N,a(t')^^CB,N,ait'))-v\,a- \ di'TTAt - t')-{cA,N A^') - ^^^B ,N A^'))) ■ (A2) 



Appendix B: Steady-state solution of the operators on the nanowire 

Expressions for the frequency functions entering the definition of the Green's function of the full system consisting 
of the nanowire and the baths, G'+(w) (cf.l^ and used in the steady state solution of the Majorana operators: 

^ , . 2i(/i-7) 2i{A + B) ij 

^lm(^) — ^ Olm -j^ 0;,m-20mod(m-|-l,4),0 0/,m-3'^mod(m,4),0 + "^<J;,m-60mod(m+l,4),0 

-r -^0Z,m-70mod(m,4),0 '';,m-20mod(m,4),0 i 'Jl,m-lOmodim+lA),Q ~r ~jr 0l,m-60niod{m,4) 



*"a a , Hf^-l) , c 2i{A^B) Z7 

^''/_,„_50niod(m-H,4),0 ~i ^ '';_,„+20mod(l+l,4),0 ^ '';,m+l0mod(l+l,4),0 ^'';,m-2''mod(m--l,4),0 

*"a a ^?!i/i^A A , 2z(A + 5) ^ i7 

^0/,m-3''mod(m-2,4),0 H <';,m+20mod(l,4),0 H 0/,m+3fmod(L4),0 ~ -^0/,m-20mod(l,4),0 

ia ^'^ X X X X ^'^ X X X X 

+ -^0/,m-l'^mod(L4),0 + ^'Ji,m+20mod(l-l,4),0 ^0/,m+l''mod(m,4),0 + ^0/,m+20mod(m,4),0 + ^'^;_„i+3dmod(m-|-l,4),0 

il X X X X ^"f X X X X ^-D^\ 

^'^i,m+6'Jmod(l+l,4),0 + -^'^i,m+5'Jmod(l+l,4),0 ^Oi,m+6<Jmod(l,4),0 ^0i,m+7'Jmod(L4),0j (t>lj 

Awn{^) = -Sj(cj)(5;^m((5;_i + (5;.2 + 5;, 3 + 61^) - ^AA'^lA^IAN-S + 5lAN-2 + ^/,47V-1 + '5;,4Ar), (B2) 
hA^) = [flL,-\{i^) ~ f]\A-^)\^niA + [vlAA - ?7L,i(^'^)]'5m,2 - i[f}L,^{i^) + ?7[^^(-w)](5™^3 - i[f]LAA + 'yl4(-'^)]'^m,4 
+ [f]R,-\{^) - j)l?„f (-'^)]'5m,4Af-3 + [tIrA^) ^ ^L^(-^^)]'5m,47V-2 ~ i[flR,-\{A + '7L^(-'^)]'^m,4Ar--l 



ilVR.lii^) + ??^J_^(-w)]<^m,4JV• 
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